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Abstract — Relative attitude formation control systems are 
developed for multiple spacecraft, based on the line-of-sight 
measurements between spacecraft in formation. The proposed 
control systems are unique in the sense that they do not require 
constructing the full attitudes of spacecraft and comparing 
them to obtain the relative attitudes indirectly. Instead, the 
control inputs are directly expressed in terms of line-of-sight 
measurements to control relative attitude formation precisely 
and efficiently. It is shown that the relative attitudes almost 
globally asymptotically track their desired relative attitudes. 
The desirable properties are illustrated by a numerical example. 

I. Introduction 

The coordinated control of multiple spacecraft in forma- 
tion has been widely studied, as there are distinct advan- 
tages [1], [2]. However, successful operation of multiple 
spacecraft requires more sophisticated control systems. For 
example, for interferometer missions like Darwin, spacecraft 
in formation should maintain a specific relative position and 
attitude configuration precisely. Relative position control and 
estimation have been well-addressed, by using Carrier-phase 
Differential GPS for precise relative navigation [3], [4]. 

Noticeable contributions on relative attitude control may 
be divided into leader-follower strategy [5], [6], behavior- 
based control [7], [8], and virtual structures [9], [10]. The 
aforementioned control systems for spacecraft attitude for- 
mation control have distinct features, but all of them are 
based on a common framework: the absolute attitude of 
each spacecraft with respect to an inertial frame is measured 
independently by using a local inertial measurement unit, 
and those measurements are transmitted to other spacecraft 
to determine relative attitudes by comparison. 

This causes restrictions on the performance of coordi- 
nated spacecraft. First, all of spacecraft should be equipped 
with possibly expensive hardware systems to determine the 
absolute attitude completely. This may increase the overall 
cost of development significantly. Second, attitude formation 
is indirectly controlled by comparing the absolute attitudes 
of multiple spacecraft in the formation. This results in a 
fundamental limitation on the accuracy of attitude forma- 
tion control systems, since measurement errors of multiple 
sensors are accumulated when determining relative attitudes. 

Vision-based sensors have been widely applied for naviga- 
tion of autonomous vehicles, where low-cost optical sensors 
are used to extract visual features to localize a vehicle [11]. 
In particular, it has been shown that line-of-sight (LOS) 
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measurements between spacecraft in formation determine 
the relative attitudes completely. An extended Kalman filter 
for relative attitude is developed based on LOS observa- 
tions [12]. The LOS measurements are also used for relative 
attitude determination of multiple vehicles [13], [14]. 

In this paper, a relative attitude formation control scheme 
is developed based on LOS measurements. Spacecraft in 
formation measure the LOS toward other spacecraft such 
that relative attitude between them asymptotically track a 
given desired relative attitude. Compared to other spacecraft 
attitude formation control systems, the proposed relative 
attitude control systems is unique in the sense that control 
inputs are directly expressed in terms of LOS measurements, 
and it does not require determining the full absolute attitude 
of spacecraft in formation or the full relative attitude between 
them. Therefore, relative attitudes are directly controlled, 
while utilizing the desirable features of vision-based sensors, 
which have higher accuracies at a relatively low cost, and 
they also have long-term stability, requiring no corrections 
in measurements as opposed to gyros. 

Compared with the preliminary work for relative attitude 
stabilization between two spacecraft [15], the control system 
proposed in this paper requires extensive analyses to take 
into full consideration of stability of time-varying systems 
for tracking, and the network structures between multiple 
spacecraft. The paper also provides stronger exponential 
stability, and numerical simulations with image processing. 

Another distinct feature of the proposed relative attitude 
control system is that it is constructed on the special or- 
thogonal group, SO(3). Attitude control systems developed 
on minimal representations, such as Euler-angles, have sin- 
gularities, and therefore their performance for large angle 
rotational maneuvers is severely limited. Quaternions do 
not have singularities, the ambiguity in representing attitude 
should be carefully resolved. By following geometric control 
approaches [16], [17], the proposed control system is devel- 
oped in a coordinate-free fashion, and it does not have any 
singularity or ambiguity. 

II. Problem Formulation 

A. Spacecraft Attitude Formation Configuration 

Consider an arbitrary number n of spacecraft in formation. 
Each spacecraft is considered as a rigid body, and an inertial 
reference frame and body-fixed frames are defined. The 
attitude of each spacecraft is the orientation of its body-fixed 
frame with respect to the inertial reference frame, and it is 
represented by a rotation matrix in the special orthogonal 



group, namely 

S0(3) = {R G 



p 3x3 



\R T R = I, detR=l}. 



Each spacecraft measures the LOS from itself toward the 
other assigned spacecraft. A LOS observation is represented 
by a unit vector in the two-sphere, defined as 



{se 



!}• 



For i,j € {1,..., n} and i ^ j, define 

Ri G SO (3) the absolute attitude for the i-th space- 
craft, representing the linear transforma- 
tion from the i-th body-fixed frame to 
the inertial reference frame, 
G S 2 the unit vector toward the j-th spacecraft 

from the i-th spacecraft, represented in 
the inertial frame, 

bij G S 2 the LOS direction observed from the i-th 

spacecraft to the j-th spacecraft, repre- 
sented in the i-th body fixed frame, 

Qij G SO(3) the relative attitude of the i-th spacecraft 
with respect to the j-th spacecraft, 

Qfj G SO(3) the desired relative attitude for Qij. 

According to these definitions, the directions of the relative 
positions in the inertial reference frame are related to the 
LOS observation bij in the i-th body-fixed frame as follows: 



Ribij j bij Ri 



(1) 



In short, bij represents the LOS observation of s^, observed 
from the i-th body. The relative attitude is given by 



Qij Rj Ri 



(2) 



which represents the linear transformation of the representa- 
tion of a vector from the i-th body fixed frame to the j-th 
body-fixed frame. Note that Qij = Qj.^. 

To assign a set of LOS that should be measured for 
each spacecraft, a graph (A^, £) is defined as follows. Each 
spacecraft is considered as a node, and the set of nodes is 
given by Af = {1, . . . , n}. The set of edges £ C Af x Af 
is defined such that the relative attitude between the i-th 
spacecraft and the j-th spacecraft is directly controlled if 
(i,j) G £■ It is undirected, i.e., (i,j) G £ (j, i) G £■ 
For each pair of two spacecraft in the edge set, another third 
spacecraft is assigned by the assignment map p : £ — > Af. As 
the edge set is undirected, the assignment map is symmetric, 
i.e., p(i, j) = p(j,i). 

For convenience, the edge set and the image of the 
assignment map are combined to form the assignment set: 

A = {(i,j,k) G £ x Af\ G £,k = p(i,j)}. (3) 

Let the measurement set Ci be the set of LOS measured 
from the i-th spacecraft, and let the communication set Cij 
be the LOS transferred from the i-th spacecraft to the j-th 
spacecraft. 

Assumption 1: The configuration of the relative positions 
is fixed, i.e., — for all i, j E Af with i ^ j. 
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Spacecraft 2 




Spacecraft 1 



Fig. 1. Formation of four spacecraft: the direction along the relative position 
of the i-th body from the j-th body is denoted by Sij in the inertial reference 
frame. The LOS observation of s;j with respect the i-th body fixed frame, 
namely bij is obtained from |TJ. 



Assumption 2: The third spacecraft assigned to each edge 
does not lie on the line joining two spacecraft connected by 
the edge, i.e., x Sjk 7^ for every (i, j, k) G A. 

Assumption 3: The measurement set of the i-th spacecraft 

is 



A = {bij,b lk £ S 2 I (i,j,fc) G A}. 



(4) 



Assumption 4: The communication set from the i-th 
spacecraft to the j-th spacecraft is given by 



C — 



{hj,bip( it j)} if (i,j)e£, 



(5) 

otherwise. 

Assumption 5: In the edge set, spacecraft are paired seri- 
ally by daisy-chaining. 

The first assumption reflects the fact that this paper does 
not consider the translational dynamics of spacecraft, and 
we focus on the rotational attitude dynamics only. The 
proposed control input does not depend on the values of 
s^, but its stability analyses is based on the first assumption 
saying that is fixed. The second assumption is required 
to determine the relative attitude between two spacecraft 
paired in the edge set from the assigned LOS measurements. 
The third assumption states that each spacecraft measures 
the LOS toward the paired spacecraft in the edge set, and 
the LOS toward the third spacecraft assigned to each pair 
by the assignment map. The fourth assumption implies that 
a spacecraft communicate only with the spacecraft paired 
with itself. The last assumption is made to simplify stability 
analysis, and the proposed relative attitude formation control 
system can be extended for other network topologies. 

An example for formation of four spacecraft satisfying 
these assumptions are illustrated at Figure [T] where 

A = {(1, 2,3), (2, 1,3), (2, 3, 1), (3, 2, 1), (3,4, 2), (4,3, 2)}. 

The measurement sets and the communication sets can be 
determined by Q and |5]l from A. For example, for the third 
spacecraft, we have £3 = {b 31 , b 32 , b 34 } , C 32 = {b 32 ,b 31 }, 
and C 34 = {b 34 ,b 32 }. 



B. Spacecraft Attitude Dynamics 

The equations of motion for the attitude dynamics of each 
spacecraft are given by 

Jitli + fij X Jiflt = m, (6) 



A. Kinematics of Relative Attitudes and Lines-of-Sight 

For any i,j 6 N, the time-derivative of the relative 
attitude is given, from (|7J, by 



(7) 



(18) 



where J; e 
and 51i g I 



p3x3 



is the inertia matrix of the z-th spacecraft, where the relative angular velocity E R 3 of the i-th 



and Ui E R are the angular velocity and 



the control moment of the i-th spacecraft, represented with 
respect to its body-fixed frame, respectively. 

The hat map A : K 3 —> so (3) transforms a vector in R 3 
to a 3 x 3 skew-symmetric matrix such that xy = x x y for 
any x, y E R 3 . The inverse of the hat map is denoted by the 
vee map V : so(3) — » R 3 . Few properties of the hat map are 
summarized as follows: 



yx -xy , 

Tf a aT\V 



x x y — xy — yx = 

tr[xA] = -tr[£(A - A T )} = -x 1 (A - A 1 ')" 

RxR T = (Rx) A , 
x ■ yz = y ■ zx — z ■ xy, 
xyz = x x (y x z) = (a; • z)y — (x ■ y)z, 
xyz — zyx — yxz, 



(8) 
(9) 

(10) 

(11) 
(12) 

(13) 



for any x,y,z e K 3 , A E M 3x3 , and R E SO(3). Throughout 
this paper, the 2-norm of a matrix A is denoted by \\A\\, and 
the dot product of two vectors is denoted by x ■ y — x T y. 
The maximum eigenvalue and the minimum eigenvalue of 
Ji are denoted by Ajv^ and A m< , respectively. 

III. Relative Attitude Tracking Between Two 
Spacecraft 

We first consider a simpler case of controlling the relative 
attitude between two spacecraft. Based on the results of 
this section, relative attitude formation control systems are 
developed later. As a concrete example, we develop a control 
system for the relative attitude between Spacecraft 1 and 
Spacecraft 2, namely Q12 = RjRi illustrated at Figure [l] 
The corresponding edge set, assignment set and measurement 
sets used in this section are given by 

£ = {(1,2), (2,1)}, 4 = {(1,2,3), (2, 1,3)}, (14) 
A = C12 = {612, b 13 }, £ 2 = C 2 i = {b 2 i, b 23 }- (15) 

Suppose that a desired relative attitude Qf 2 (t) is given as 
a smooth function of time. It satisfies the kinematic equation: 



Ql2 — Ql2^12 



d r\d 



(16) 



where £lf 2 is the desired relative angular velocity. Note that 
these also yield Q21 = (Qi2) T from ([2J, and it satisfies 



?21 — Q21^2lJ 



(17) 



where = -Qf 2 0f 2 . 

The goal is to design control inputs u\ , U2 in terms of the 
LOS measurements in C\ \JJH2 such that Q12 asymptotically 
follows Qf 2 , i.e., Quit) —> Qf 2 (t) as t -t 00. 



spacecraft with respect to the j-th spacecraft is defined as 

Qij = fii - Qjjflj. (19) 

From ([TJ and (|7j, the time-derivative of the LOS measure- 
ment bij is given by 



bij Ri Sij QiRi &ij — bij X 



(20) 



Let bijk € K be bijk = bij x b^. From (20 1 and (12i, it 
can be shown that 

bijk = (bij x Qi) x bik + hj x (h lk x 

= —(fti ■ b tk )b t .j + (CI, ■ bij)b ik = bijk x Vti. (21) 

B. Relative Attitude Tracking 

It has been shown that four LOS measurements in L\ U 
£2 completely determine the relative attitude Q 12 from the 
following constraints [15]: 

bl2 = — Ql2^21; 
^123 <2l2^213 



\bi 



(22) 
(23) 



23 



'213 



These are derived from the fact that four unit vectors, namely 
{S12, S13, S21, S23} lie on the sides of a triangle composed 



of three spacecraft. The first constraint (22i states that the 
unit vector from Spacecraft 1 to Spacecraft 2 is exactly 
opposite to the unit vector from Spacecraft 2 to Spacecraft 
1, i.e., S12 = —821. The second constraint (23 1 implies 
that the plane spanned by s 12 and s 13 should be co-planar 
with the plane spanned by S21 and S23. These geometric 
constraints are simply expressed with respect to the first 
body-fixed frame to obtain ( f22"l > and ( f2"3"] >. The relative attitude 
Q12 is determined uniquely by the LOS measurements 
{bi2, b\3, b 2 i, 623} according to (22 1 and (23 1. 



We develop a relative attitude tracking control system 
based on these two constraints. More explicitly, control 
inputs are chosen such that two constraints are satisfied when 
the relative attitude is equal to its desired value. As both 
constraints are conditions on unit vectors, controller design 
similar to tracking control on the two-sphere. From now on, 



variables related to the first constraint (22 1 (resp., the second 



constraint ( 23 1) are denoted by the sub- or super-script a 
(resp., P). 

First, configuration error functions are defined as 

\\\b 2 i + Q? 2 M 2 = 1 + 6ai ■ Q? 2 &i2, (24) 

1 + 



*?2 



12 



a\2 



'213 



Qt 2 b 



12°123, 



(25) 



where a\2 = »2i 



I&213IIH&123II e M. Since ||b 



\hj x b ik \ 



\Rjsa x Rjsi 



'ijk\\ — 

\sn x Sjfcll, the constant 



a\2 is fixed according to Assumption [T] and it is non-zero 
from Assumption [2] Next, we define the configuration error 
vectors as 



e? 2 = (Q&&21) x 612, 



<%i = (Qia&ia) X 621, (26) 



ai2 



-(<&&2is) x 6 



123j 



«21 



(Q126123) X 6 2 



13- 



(27) 



As 612,621 are unit vectors, and from the definition of 

ai2, a 2i> we can show that ||e" 2 ||) II e 



21 » ll e 12ll» I' 21 



We also define the angular velocity errors: 



e ai = n t - n?, 



en 2 



— r^n 



< 1. 



(28) 



where the desired absolute angular velocities QJ, are 
chosen such that 



n? 2 (*) = n?(t)-Q^(t)n^(t). 



(29) 



Any desired absolute angular velocities satisfying (29 1 can 
be chosen. For example, they can be selected as 

Using these desired angular velocities, the derivative of the 
desired relative attitude can be rewritten as 



Ql2 ~ Ql2^1 ^2Ql2- 



(30) 



It is assumed that the desired angular velocities are bounded 
by known constants. 

Assumption 6: For known positive constants Bd, 

||fi?(t)|| <B d , \\Q*(t)\\ < B*, 

for all t > 0. 

The properties of these error variables are summarized as 
follows. 

Proposition 1: For positive constants k" 2 ^ fcf 2 , define 

*i3 = fcia*ia + (3D 

(32) 

(33) 



C12 — «-12 c 12 1 /t 12 e 12i 



£21 — ^21 e 21 + ^21 e 21i 

where fc 21 = A; 21 , fc 21 = fcf 2 . The following properties hold: 

(i) ei2 = -Qi x e 2 i, and ||ei 2 || = ||e 2 i||. 
(") S*i2 = e i2 ' e Oi + e 2 i • en 2 . 

U/3 



(iii) ||ei 2 || < (fcf 2 + /c^)(||eo 1 1| + ||en 2 
||e 2 i||<(A:f 2 + fcf 2 )(||ea 1 || + ||eo 2 



S d ||ei2||, 
S d ||e 2 i||. 

(iv) If ^i 2 < ip < 2min{fc" 2 , fcf 2 } for a constant ip, then 
is quadratic with respect to ||ei 2 ||, i.e., the following 
inequality is satisfied: 

V\Jei2|| 2 <*12 <^i 2 ||ei2|| 2 , (34) 



4\ 



where the constants tp 12 ,ipi 2 are given by 

= min{fc^ 2 ,fcf 2 } 

2max{(fcf 2 )2, (fcf 2 - fcf 2 )2} + 2(fc« 2 + ' 

min{fcf 2 ,fef 2 }(fef 2 + fcf 2 ) 



i{(fcf2) 2 ; (fcf2) 2 }(2min{^ 2 ,fcf 2 }-V)' 



Proof: From ( |2"6*j i, ( [TO} , e" 2 is given by 

e? 2 = (Qiihi) x 612 = Q^i{6 2 i x Q? 2 6i 2 } = -Q&ejj. 

Similarly, ef 2 = —Qii e 2i- Th ese show (i). 

Since a; T T/ = tr[a;j/ T ] for any x, y £ M 3 and from (Jlj, the 
configuration error function can be written as 



* 12 = K 2 + k 



^12, 



12 t- fc?2tr[6 2 i6( 2 Qf 1 ] + ^tr[6 2 i 3 "i23 



-^123 ( 92l] 

Ctl2 

tr[i^(fef 2 « 2 i«f 2 + — 52133^3)^^21] 

Ol2 

(35) 



— k" 2 + k 1 



^k^ + k^-t^K^Q^], 



where s ljk £ 
Since s i2 = -S21, 
rewritten as 



denotes s 

S213 



sa x s ife , and K 12 £ 



p3x3 



^J^y, the matrix Xi 2 can be 

S123 



ifl2 = fc" 2 Si 2 sf 2 



ft 12 



Sl23| 



T 

2 S 123S 123 . 



(36) 



The derivative of the configuration error function with 
respect to R\ along the direction of SRi = Rifji for a vector 
r/i £ M 3 is given by 



D_R 1 *i2(-Ri, -R2, Q12) ■ R1V1 



de 



4-i 2 (i?iexpe?7i,i?2,Qi 2 ) 



e=0 



= -tr[R^K 12 RimQ d 21 } = -tvlfnQ^RlKiiRx] . 
Using Q and ([8j, this is rewritten as 
U Rl ^ 12 (R 1 ,R 2 ,Q d 12 ) ■ Rxm 

= Vi ■ [Q d 21 RlK 12 R x - RjK 12 R 2 Q( 2 y (37) 

k 

= »ft • (-fcf 2 <9li^2 «ia x R[s 12 - -^Q d 2X R T 2 s 12Z x Rjs 123 ) 

a X2 

= Vl ■ (fc? 2 Q21&21 X 6i 2 + ^g^6 2 i 3 X i?T &123 ) 

Ol2 



Vi ■ e 12 . 



(38) 



In short, the error vector e\ 2 is the left-trivialized derivative 
of with respect to R x . Similarly, we can show that 

D J R 2 *12(i?l,i?2,<3^) ' ^2??2 = m ■ e 2U 

U Q ^ 12 (R 1 ,R 2 ,Qf 2 )-Qf 2 fi 12 = -vi2-e 12 . 

Therefore, the time-derivative of the configuration error 
function is given by 



^12 = ei2 • ill + e2i • ^2 - ei2 ■ ft v 



12- 



Substituting (29 1 and since Qf 2 e i2 = —e 2 \, we obtain 

*i2 = cia • - nf) + e 21 ■ (o 2 - ng), 

which shows (ii). 

Next we show (iii). From ( [37] ), ei 2 can be written as 

ei2 = {Q&xBZKhRx - R[K 12 R 2 Qf 2 ) v . 



(39) 



From attitude kinematic equations (j7]), ( [30| , the time deriva- 
tive of the error vector ei2 is given by 

ei2 - {Q^RjKuRiflt + ^RjK^R^) 



- Q&tMI tfiafliQai + Q d i2R T iK 12 R 2 fl 2 )Qi 2 

- {ClfQ^RlKuRx + RjK 12 R 2 Qi 2 nf) 

+ Q 21 (tl 2 R 2 K 12 RiQ 21 + Q ' 2 J Rj K \ 2 R 2 tl 2 )Qf 2 . 

Using ( f2"8j ), this is rewritten as 

I12 = (Q^RlK^Riin, + e ni R[K 12 R 2 Q d 12 ) 
+ (Q 21 R 2 K 12 RiClf + niR]K 12 R 2 Qi 2 ) 

- Qii{^n 2 RiK\ 2 RxQ 2l + Q n R 1 K 12 R 2 ea 2 )Q 12 

- {nfQ d jR 2 ~K 12 R 1 + R~[K 12 R 2 Q d l2 ni). 



Substituting ((36b, and applying ((8j), (13 1, this reduces to 



era = K 2 {-Qi x b 21 x (e 0l 6 12 ) + &12 x Q2i(eo 2 &2i)} 

+ — {-Q21 fo 213 X (e 0l &123) + &123 x Qfj (en 2 6213)} 
^12 

-Afeia, (40) 

which shows the first inequality of (iii) from Assumption [6] 
The second inequality of (iii) can be shown similarly. 

Next, to show (iv), we use the following properties given 
in [15]. For non-negative constants fi,f 2 ,f 3 , let F — 
diag[/i, f 2 ,f 3 ] G K 3x3 , and let P G SO(3). Define 

$ = -K[F(I - P)], (41) 



e P = i(FP-P T F) v , 



(42) 



Then, $ is bounded by the square of the norm of ep as 



kl >\e P \\ 2 <<P< 



ep 



(43) 



h 2 + /13 "~ J " ' hsQii - 4>) 
if $ < < hi for a constant 0, where /i, are given by 

hx = min{/i + / 2 , ,/ 2 + / 3 , / 3 + A}, 
h 2 =tmx{(f 1 - h)\ {h-h)\ (iWi) 2 }, 
/ l3 =max{(/ 1 + / 2 ) 2 , (/ 2 + / 3 ) 2 , (/3 + /1) 2 }, 
/i 4 = max{/i + f 2 , f 2 + f 3 , f 3 + fx}, 
/i 5 =min{(A + / 2 ) 2 , (/ 2 + / 3 ) 2 , (/ 3 + /i) 2 }- 
We apply this property by showing that the configuration 



error function vj/ 12 given at (35 1 can be rewritten as (41 



At (36 1, the matrix K\ 2 can be decomposed into Ki 2 = 
UGu* , where G is the diagonal matrix given by G = 
diag[fef 2) fcf 2 ,0] G M 3 , and U is an orthonormal matrix 



defined as U = [si 2 , 



G SO(3). Using the 



||si23ll ' UsiaXsiasll 

property, tr[AB] = tr[BA] = tr[B T A T ] = tr[A T B T ], the 
configuration error function can be rearranged as 

*ia = tr[G(/ - U T RxQ^RlU)}. (44) 

Therefore, if we choose F = 2G and P = XJ T R x Q d 21 R\\J , 



we obtain = <&• Also substituting these into (42 1, 

e P = GV T R 1 Q d 21 RlV - U T R 2 Q d 12 R^UG 
= U 7 \K X2 R x Q d 21 Rl - R 2 Q d 12 RjK 12 )U 
= U T R2Q12 (Q 2X R 2 K\ 2 Ri — R x Ki 2 R 2 Q 12 )Q 21 R 2 U. 



From this, we have e P = U T R 2 Q d 2 e 12 by (39i, (lOi 
Therefore, ||ep|| = 1 1 e 12 1 1 . Then, (43 1 yields (iv) with fx 
2fc? 2 , f 2 = 2fcf 2 , f 3 = 0. 



Using these properties, we develop a control system to 
track the given desired relative attitude as follows. 

Proposition 2: Consider the attitude dynamics of space- 
craft given by ([6j, (7jl for i G {1,2}, with the LOS measure- 
ments specified at (14). A desired relative attitude trajectory 



is given by (I61. For positive constants kf 2 ^ k^ 2 ,k 21 



kf 2 , k 21 — ki 2 , kfi 1 , k(i 2 , control inputs are chosen as 



ka i en i + ClfJi(en i + flf) + JO 



(45) 



where G £. Then, the following properties hold: 

(i) There are four types of equilibrium, given by 
the desired equilibrium (Q,£li 2 ) = (Qi 2 ,^t 2 ), 
and the relative configurations represented by 
Q d 2 = R 2 r UDU T R 1 and Q 12 = ftf 2 where D G 
{diag[l, -1, -1], diag[-l, 1, -1], diag[-l, -1, 1]} 
and U G SO(3) is the matrix composed of eigenvectors 
of K12 given at ( [36| . 

(ii) The desired equilibrium is almost globally exponen- 
tially stable, and a (conservative) estimate to the region 
of attraction is given by 



E 

z=l,2 



ttia(0) <^<2min{fc«,fcf 2 }, (46) 
Am, || eo^O)!! 2 <2(V>-*i 2 (0)), (47) 



where ip is a positive constant satisfying ^ < 
2 min{fc" 2 , fcf 2 }, and Am ; denotes the maximum eigen- 
value of Ji. 

(iii) The undesired equilibria are unstable. 

Proof: From ([6}, ( |45] >, and rearranging, the time- 
derivative of Jieo i is given by 



Jiea t = -(eni + Of) x J^en, 



■ 0?) + Uj - 



(48) 
(49) 



The equilibrium configurations are where ei 2 = e 2 i = 
eui = = 0, which corresponds to the critical points 



of the configuration error function given by (44i. In [17], it 
has been shown that there are four critical points: 

RiQtxRl G {I, UDxU T , UD 2 U T , UD 3 U T }, 

where D x = diag[l, -1, -1], D 2 = diag[-l, 1, -1], D 3 = 
diagj— 1, — 1, 1]. This shows (i). 

Next, we show exponential stability of the desired equi- 
librium. A sufficient condition on the initial conditions to 



satisfy (34i is obtained from the following variable: 



1 



en 2 ' J2£n 2 



From ( |49] > and the property (ii) of Proposition [T] the time- 
derivative of U is simply given by 



W = -fcn 1 ||en 1 || 2 -fco 2 ||eo 2 || 2 <0, 



(50) 



which implies that U(t) is non-increasing. For the initial 
conditions satisfying ( |46| , ( |47] i, we have 

W(0) < - E A M J|en s (0)|| 2 + *i 2 (0)<^. 

i=l,2 

As (0) is non-increasing, 

*i 2 (t) < W(t) < 1/(0) <4>< 2min{^ 2 ,fcf 2 }. 

Therefore, ^12 (i) <ip<hi for all f > 0, and the inequality 
( |34| i is satisfied. 

Let a Lyapunov function be 

V = U + c(ei2 • + e 2 i • en 2 ), 
for a constant c. Using ( [34| , it can be shown that 

E ^iMijZii ■ V- ^ z^M ijZij , (51) 

,||enj|] € R 2 , and the matrices 



where zy = [||e 
Mij,Mij e M 2x2 are defined as 



M, 



— y 

-c A r 



Mi 



c A 



(52) 



for (i,j) € £ = 1(1, 2), (2, 1)}, where it is assumed that 

ip. . — ip .., — From (50 1, we obtain 

— y — ~* 



V= ^ -fej|enj| 2 + c(e„- • en, + e i:) - • eoj. (53) 



From (49 1, and using the fact that He^l < kfj + k^, we have 



en s • ey < --L||ey|| 2 + J* (fc» + fc^, , «, 



■J3- 



Am< A r , 
^M t B d + kn .. 

A l|£ 



(54) 



Together with the property (iii) of Proposition [T] this yields 
the following inequality of V: 



v= E -fa 



A 



c (fc«+fc£)(l + ^))|| ent || a 



ijes 
c 

A M; 

c 



lleyf + c(fcg + *£)l|enJ||e n J 
((A Mi +A m4 )S d + fe i )||e«||||enJ|. 



(55) 



This can be rewritten as the following matrix form: 

V < -zf 2 W 12 z 12 - z 2 r 1 W 2 iz 2 i - CS^12Cl2, 

where £x 2 = [Il e ni || , ||en 2 ll] T ^ M 2 and the matrices 
W12, W21, Y 12 <= M 2x2 are given by 

xfr -^(A^ + fcaJ- 



(56) 



lis 



- Cfcl 2 (l + ^) 

-2cfci2 fcn 2 



-2cfc 



12 



cfc 12 (l+^) 



(57) 



for £ £■ Here, A, denotes A^ 



denotes k 



y 



fc£. It can be shown that if the constant c 



Am s + A J( 



and k 



1 j 



is sufficiently small, then all of the matrices My, Wi ; 



for (i, j) € £, and Y\ 2 at (51 1 and (55 1 are positive definite. 
For example, My is positive definite if c < \J^^Knv 
This shows that V is positive definite and decrescent, and 
V is negative definite. Therefore, the desired equilibrium is 
exponentially stable. 

Next, we show (iii). At the first type of undesired equilibria 
given by RiQ^R^ = UDiU T and e Ql = en 2 = 0, the 

v i2- 



value of the Lyapunov function becomes V = 2fcf - Define 



W = 2fcf 2 - V. 
Then, W = at the undesired equilibrium, and we have 



^||eo I || 2 + C ||e l3 ||||eo 1 ||} + (2fcf 2 -vI/)<W. 



- Ei 

Due to the continuity of ^?, we can choose i?i and R 2 
arbitrary close to the undesired equilibrium such that (2fc 2 — 
if?) > 0. Therefore, if ||enj| is sufficiently small, we obtain 
W > 0. Therefore, at any arbitrarily small neighborhood 
of the undesired equilibrium, there exists a set in which 
W > 0, and we have VV = -V > from (55 1. Therefore, 
the undesired equilibrium is unstable [18, Theorem 3.3]. 
The instability of other types of equilibrium can be shown 
similarly. This shows (iii). 

The region of attraction to the desired equilibrium ex- 
cludes the union of stable manifolds to the unstable equi- 
libria. But, the union of stable manifolds has less dimension 
than the tangent bundle of the configuration manifold. There- 
fore, the measure of the stable manifolds to the unstable 
equilibria is zero. This implies the desired equilibrium is 
almost globally exponentially stable [16], which shows (ii). 

■ 

This states that almost all solutions of the proposed control 
system, excluding a class of solutions starting from a specific 
set that has a zero-measure, asymptotically track given the 
desired relative attitude. As the control inputs are expressed 
in terms of observations, in addition to angular velocities, and 
the full relative attitude does not have to be constructed at 
each time. These results can be considered as a generalization 
of the preliminary work in [15], but it is a nontrivial 
extension as the several properties of the error variables 
should be considered to show a stronger exponential stability 
for tracking problems. 

IV. Relative Attitude Formation Tracking 

The relative attitude control system between two space- 
craft developed in the previous section can be used as a 
building block for a relative attitude formation control system 
for multiple spacecraft. In this section, we generalize it for 
daisy-chained relative attitude formation control network. 

A. Relative Attitude Tracking Between Three Spacecraft 

We first consider relative attitude formation tracking be- 
tween three spacecraft, given by Spacecraft 1, 2, and 3, 



illustrated at Figure [T] The corresponding edge set and the 
assignment set used in this subsection are given by 

£ = {(1,2), (2,1), (2, 3), (3, 2)}, (58) 
.4 = {(1,2, 3), (2, 1,3), (2, 3,1), (3, 2,1)}. (59) 

For given relative attitude commands, Qf 2 (t) 7 Q 23 (t), tne 
goal is to design control inputs such that Q\ 2 (t) — > Qf 2 (t) 
and Q 23 (t) ->• Q d 3 (t) as t — > oo. 

The definition of error variables and their properties de- 
veloped in the previous section for two spacecraft are readily 
generalized to any (i,j, k) £ A in this section. For example, 
the kinematic equation for the desired relative attitude Q 23 
is obtained from ( fT6| > as 

Q23 = Q23^23> 

where il 23 is the desired relative angular velocity. Other 
configuration error functions and error vectors between 
Spacecraft 2 and Spacecraft 3 are defined similarly. 

The desired absolute angular velocities for each spacecraft, 
namely fif , f2 2 , and il 3 should be properly defined. For the 
given ilf 2 , f2 23 , they can be arbitrarily chosen such that 

n d 12 (t) = n d 1 (t)-Qi 1 (t)n d 2 (t), (60) 

n&(t) = ng(«) - QUmi(t). (6i) 

For example, they can be chosen as flf = £lf 2 , f2 2 = 0, 
^3 — ~Q 23 ^ 23 - Assumption |6] is considered to be satisfied 
such that each of the desired angular velocity is bounded by 
a known constant B d . 

Proposition 3: Consider the attitude dynamics of space- 
craft given by d6j, ^ for i G {1,2,3}, with the LOS 
measurements specified at (|59|. Desired relative attitudes are 



given by Qf 2 (t), Qi^t). For positive constants k" j ,k^ j ,kQ i 
with fc« ± fc£, fcg = kf t , fcg = fcf. for (i, j) e f , 

ui = -e 12 - fc fil e fil + n? J 1 (e ni + fif) + J(%, (62) 

u 2 = -^( e 2i + e 23 ) - fen 2 en 2 + ft!} J 2 (en 2 + ng) + JQ£, 

(63) 

u 3 = -e 32 - ka 3 e Q3 + tlp 3 (eQ 3 + O3) + ^3> (64) 

Then, the desired relative attitude configuration is almost 
globally exponentially stable, and a (conservative) estimate 
to the region of attraction is given by 

*i 2 (0) + * 23 (0) < V < 2min{A&, (65) 
A Ml IK (0)|| 2 + 2A M2 ||en 2 (0)|| 2 + A Ma |)en 3 (0)|| 2 
< 2(ip - ¥ 13 (0) - #aa(0)), 



(66) 



where -0 is a positive constant satisfying ip < 
2min{fcf 2 , k 12 , k 23 , k 23 }. 

Proof: The time-derivative of J\C\ and J 3 e 3 are given 

^2 e 2 is given by 

(e 2i + e 23 ) - fcfi 2 en 2 - 



by (49 1, and the time-derivative of J 2 e2 is given by 

1 
2 



hen 



(J 2 en 2 + J^i)en 2 



(67) 



Define 



1 



en 2 ■ Ji&o. 2 + -en 3 ' J3 e n 3 



+ *i2 + # 23 . 

From ( |49] i, |67|, we have 

U = -fenjenj 2 - 2fcn 3 l|efi 2 || 2 - fen 3 ||en 3 || 2 , 

which implies that U(t) is non-increasing. For the initial 
conditions satisfying |65] l and (661, we have W(0) < tp. 
Therefore, 

#12(0 + # 23 (<) < W(t) < W(0) < V < 2 min{fc? 2 , fcf 2 , fc 23 , / 



Therefore, the inequality ( |34] i holds for both of #i 2 and # 2 3. 
Let a Lyapunov function be 

V = U + cen± ■ ei2 + ce^ 2 ■ (e 2 i + e 23 ) + ce fi3 ■ e 32 . (68) 



From ( 34 1, we can show that this Lyapunov function satisfies 



the inequality given by ( 5 1 1. The time-derivative of the 
Lyapunov function is given by ( |53~| . 

For £ {(1, 2), (3, 2)}, the upper bound of e^. • eu is 
given by p4) i. From (67 1, the upper bound of in 2 ■ (e 2 i +e 23 ) 
is given by 

1 „ 



e o 2 



(e 2 i + e 23 ) < ||e 2 i 

Am 2 



w 2 



A, 



(fcai + fc 2 3)||en 2 || 



62311 
^M 2 B d 



An 



e 2 i + e 23 ||||ef2 2 |[ 



The upper bounds of ||e 32 || and 1 1 e x2 1 1 are given by the 
property (iii) of Proposition [T] Additionally, using (40i, we 
can show that 

P21 +e 23 || < (A a i + A 2 3)(||en 1 || + 2||en 2 || + ||en 3 ||) 
+ B d \\e 2 i +e 23 ||. 

Applying these bounds to the expression of V and rearrang- 
ing, we obtain This can be written as a matrix form as 



V < -z[ 2 W 12 z X2 - z 213 W 213 z 213 - z 32 W 32 z 32 

— Coital C21 — C23^23C23, 



(69) 



where the matrix Wi 2 ,W 23 £ M 2x2 are given as (56i, and 

I 2 . The matrices W213, Z 2 i, 



z 2 \ 3 = 

Z 23 £ 



[\\e21 + e 23 \\,\\en 2 \\\ £ 
|2x2 are (jefjjjej as 



IT"; 



213 



Aw 



-j^(X 2 B d + kn 2 ) 
CX 2 B d + k n2 ) 2k n2 - ck 213 (2 + ^ 



-"2 



^21 — ~r 



^23 — 



A 

2fco 2 - cA- 2 i 3 (2 , A 
-c(6fci2 + 2fc 23 ) 

2fco 2 - cA: 2 i3(2 t x 
-c(6fc 23 + 2fc 12 ) 



M-2 



-c(6fc 12 + 2fc 23 ) 
fe^ - cfei 3 (l + ^ 



-c(6fe : 



23 



2fe i2 ) 
c fel3(l + 



where ^213 = ^12 + ^23- It can be shown that if the constant c 



is sufficient small, all of matrices at (51 1 and (69 1 are positive 
definite, which implies 

V<-A m (M/ 12 )(|| ei 2|| 2 + ||eo 1 || 2 ) 



X m {W 213 )(\\e 21 + e 23 \\ 2 + \\e n2 \\ 2 ) 
A m (W 32 )(||e 31 || : 



e ^3 



1 



<-A m (^ 12 )|-(||e 12 || 2 + || e21 || 2 ) + || eni || 2 
- \ m (W2iz)\\en 2 \\ 2 



A m (^ 32 )<j-(|| e32 || 2 + ||e 23 || 2 ) + ||en 3 || 



where A m (-) denotes the minimum eigenvalue of a matrix, 
and we use the fact that ||ei 2 || = ||e 21 ||, ||e 23 || = ||e 32 ||. 
Therefore, the desired equilibrium is exponentially stable. 

To show almost exponential stability, it is required that the 
fifteen types of the undesired equilibria, corresponding to the 
critical points of ^i 2 and \&23> are unstable. This is similar 
to the proof of the property (iii) of Proposition |2j and it is 
omitted. ■ 

The control inputs for Spacecraft 1 and Spacecraft 3 at the 
both ends of graph are identical to (45 1 at Proposition [2] The 
control input for Spacecraft 2, which are paired with both 



of Spacecraft 1 and 3, is also similar to (45 1 except that the 
configuration error vectors for Spacecraft 2, namely e 2 i and 
e 23 are averaged. These ideas can be generalized to relative 
attitude formation tracking between an arbitrary number of 
spacecraft as follows. 

B. Relative Attitude Formation Tracking Between n Space- 
craft 

Consider a formation of n spacecraft, i.e., Af = 
{1, . . . , n}. According to Assumption^ spacecraft are paired 
serially in the edge set. For convenience, it is assumed that 
spacecraft are numbered such that the edge set is given by 

£ = {(1,2),..., (n - 1, n), (2, 1), . . . (n, n - 1)}. (70) 

The assignment set is given by ([3]) for an arbitrary assignment 
map satisfying Assumption [2] The desired relative attitudes 
Qfj for £ £ are prescribed. The definition of error 

variables and their properties developed in Section III are 



generalized to any (i, j, k) £ A. The desired absolute angular 
velocities Clf for i £ Af are chosen such that 

n&(t) = 0*(t) - Q%m%t) for (ij) £ £. (71) 

Proposition 4: Consider the attitude dynamics of space- 
craft given by {6}, (|7| for i £ {1, ...,n}, with the LOS 
measurements specified by f70|),(|3]l. Desired relative attitudes 
are given by QfAt) for £ £. For positive constants 

fcg,fc£.,*n 4 with k% + k%, 
inputs chosen as 

ui = -ei 2 - kn 1 eQ 1 



0-1.0 



k?i, the control 



OiJi(e ni + fii) + Jflt 



-ei2 

'\( e P, P -i + e P, P +i) ~ k n p en p + (l d p J P {ea p 
- Jtl d p for p e {2, . . . , n - 1}, 
-e n ,„_i - kn n en n + (l d n J n {en n + OJl) + J(l£, 



(72) 
(73) 




Fig. 2. Relative attitude formation tracking for 7 spacecraft: the lines-of- 
sight measured by each spacecraft are denoted by arrows, and the dotted 
line between two spacecraft implies that they are paired at the edge set. 



Then, the desired relative attitude configuration is almost 
globally exponentially stable, and a (conservative) estimate 
to the region of attraction is given by 

n-l 

E *M+i(0) < V> < 2 nun {A$ +lJ 

* — 4 Kt<n — 1 ' 



A Ml ||en 1 (0)|| 2 + 2 V A M4 ||e ni (0)|| 2 + A M J|en n (0)|| 



i=2 



<2(^ 



n-l 

E 

i=l 



(0)), 



< 



where ip is a positive constant satisfying rp 
2mini< 4 <„_i{fc^ +1 , k P l i+1 }. 
The proof of this proposition is a straightforward, but tedious 
extension of the proof of Proposition|3] Due to page limit, the 
detailed proof is omitted, but numerical results for multiple 
spacecraft are provided in the next section. 

V. Numerical Example 

Consider the formation of seven spacecraft illustrated at 
Figure [2] The corresponding edge is given by ( |70"| i with n = 
7, and the assignment set is 

A = {(1, 2, 3), (2, 1, 3), (2, 3, 4), (3, 2, 4), (3, 4, 5), (4, 3, 5), 
(4, 5, 7), (5, 4, 7), (5, 6, 7), (6, 5, 7), (6, 7, 5), (7, 6, 5)}. 



The desired relative attitudes for 



Q d 3i 



and Qf 5 are 

given in terms of 3-2-1 Euler angles as Qfaif) = 
Qf 4 (a(t),/3(t), 7 (t)), Qf 5 (t) = QU<t>(t), where 

a(t) = sin0.5i, f3(t) = 0.1, j(t) = cost, 
(j)(t) = 0, 6{t) = -0.1 + cos 0.2*, V (*) = 0.5 sin 2*, 

and Qf 2 (*) = Qi 3 (t) = Qf 6 (*) = J, Qg 7 (t) = (Qf 5 (i)) T . It 
is chosen that ilf(t) = 0, and other desired absolute angular 



velocities are selected to satisfy (71 



(74) 



The initial attitudes for Spacecraft 3 and 6 are chosen 
as R 3 (0) = exp(0.9997rei) and R 6 (0) = exp(0.9907re 2 ), 
where ei = [l,0,0] T ,e 2 = [0,1, 0] T G M 3 . The initial 
attitudes for other spacecraft are chosen as the identity 



(a) Relative attitude error functions 
*12,*S3,...,*67 



0^ 



(c) Relative angular velocity error 



(b) 



Relative attitude error vectors 

2 ' e QlZ ' * * " ' e Q67 




(d) Control moments iti 
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Fig. 3. Numerical results for seven spacecraft in formation (blue, green, 
red, cyan, magenta, and black in ascending order) 



matrix. The resulting initial errors for the relative attitudes 
Q23 and Q67 are 0.997rrad = 179.82°. The initial angular 
velocity is chosen as zero for every spacecraft. 

The inertia matrix is identical, i.e., J. L — diag[3, 2, 1] kgm 2 
for all i 6 J\f. Controller gains are chosen as fc^. = 7, 
feg. = 25, and fcg = 25.1 for any <= £. 

Tracking errors for relative attitudes and control inputs are 
shown at Figure [3] where the relative attitude error vectors 
are defined as e Qz] = \{{Q%f Qij - QjjQfj)" € M 3 . These 
illustrate good convergence rates. 
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